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. Abstract 

f-n | We present an asymmetric g-Painleve equation. We will derive this using q-orthogonal 

polynomials with respect to generalized Freud weights: their recurrence coefficients will 
obey this g-Painleve equation (up to a simple transformation). We will show a stable 
method of computing a special solution which gives the recurrence coefficients. We estab- 
lish a connection between the newfound equation and a-q-Py. 

> 

. 1 Introduction 

On 
O 

The aim of this work is to identify new g-discrete Painleve equations. What exactly are the 
discrete Painleve equations (d-P)? One could state that a d-P is a second-order, nonautonomous 

OO ' 
O 



integrable mapping having one of the celebrated Painleve equations for a continuous limit. This 
description does not quite capture the whole story, so let's look at some important notions about 
discrete Painleve equations. 
Although some discrete Painleve equations were found avant la lettre, the first steps were made 
by Brezin and Kazakov pQ who found an equation (now known as d-Pi) and computed its 
continuous limit, the continuous Painleve I. The definition of singularity confinement [5] was 
maybe the most important step in the evolution of the field: it is the discrete analogue of the 
Painleve property and describes the behaviour of singularities throughout the evolution of a 
discrete equation: 

Definition 1.1 (Singularity Confinement Property). Consider a difference equation 
with independent variable n and dependent variable x n . If x n is such that it gives rise to a sin- 
gularity for x n+ i, then there exists ap EN such that the singularity is confined to x n+ i, . . . , x n+p 
and x n+p+ i depends only on x n -i,x n -2, ■ ■ ■ 

Using only this criterion of singularity confinement, discrete analogues for (continuous) Pm, Piv 
and Py were identified [5], [9]. It's striking that in the analogue of Pm we see the independent 
variable n entering the equation only in an exponential way: 

(x n + a)(x n + (3) 



■^n+l^n—l 



(jq n + l)(Sq n + 1) 
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with a, j3, 7, 5 constants. Therefore, it is considered to be a so-called q-discrete equation as the 
nature of the equation is multiplicative rather than additive like, e.g., d-Pi 



z w + 7 (-l)" ^ x 

■En— l ~T X n -+- X n +i — h 

Xn 

with z n = an+(3 and a, /3, 7, 5 constants. This equation shows an odd-even dependence through 
the factor (—1)™ when 7 7^ 0. In this case we can introduce new variables u n = x 2n , v n = x 2n +i 
that lead to the system 

Vn _ 1 + Un + Vn = ^±M±^ + 5 = *-^ + 5 
u n + v n + u n+1 = Mnj^zl + s = ^2 + 6 

with z n = 2na + f3. This system is known as a— d-Pi, an asymmetric discrete Painleve equation. 
A list of a few important discrete Painleve equations was compiled by Peter Clarkson, and can 
be found in [IT] . See also [I] for an overview of discrete Painleve equations. 



The link between orthonormal polynomials and discrete Painleve equations is well established 
(0) [Z], [II])- Given a positive measure w on a set A C R and assuming all the moments for 
w exist, i.e. | j A x k w(x)dx\ < 00, we denote with {p n } the set of orthonormal polynomials with 
respect to w: 

p n {x)p m {x)w(x)dx = 5 mn m,n>0. 
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These orthonormal polynomials are unique if we choose the leading coefficient to be positive. 
Orthonormal polynomials satisfy a three-term recurrence relation of the form 

xp n (x) = a n+ xp n+1 (x) + b n p n (x) + a n p n „i(x) n>0 

with p_i = 0. Here we can find the recurrence coefficients a n and b n as the coefficients showing 
up in the Fourier expansion of xp n (x): 

Xp n (x)p n _i(x)w(x)dx, (1.1) 

xp n {x)w{x)dx. (1.2) 



.4 

If we write p n {x) = •y n x n + . . ., we can find a n = by comparing the leading coefficients of 
both sides of the recurrence relation. Since we chose to take positive leading coefficients, a n 
will be positive as well. 

During the last few decades, many generalizations of well-known weight functions have been 
shown to give rise to discrete Painleve equations for the recurrence coefficients a n and b n . 
As far as we can see, the very first appearance of a d-P was in the context of orthogonal poly- 
nomials, when Shohat [10J found a nonlinear recurrence relation for the recurrence coefficients, 
which is now known as d-Pi. Fokas, Its and Kitaev found the connection of d-Pi to Freud 
weights in [2j. Nijhoff [8] found a g-discrete Painleve equation in the context of orthogonal 
polynomials: he considered a g-generalization of the Hermite polynomials on the exponential 
lattice and found non-linear recurrences of order >2. A second order recurrence relation was 
found by the third author [TT] which he called q-Pi- We will study a slight extension of this 
Freud weight to recover asymmetric g-Painleve equations. 
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In particular we consider in Section 3 the weight 

{q 4 x 4 ;q 4 ) OQ \x\ a 

w(x) = — ; tt — T . — a > —I 

K J (1 - q A ) a / 4 

on the g-exponential lattice L = {±g n , n G N}, with q G (0, 1). We will show that b n = and 
Vn — a nq 1 ~ n satisfies 



q n a {y n yn+i + q a )(y n y n -i + q c 



q a — q a y 2 n even 
1 — y 2 n odd 



which we will call q-P/. The case a = was already obtained in [TT]. In Section 4 we consider 
the more general weight 

\x\ a {q 2 x 2 ]q 2 ) OD {cq 2 x 2 ]q 2 ) 00 

W(X) = ; — -r. , a > — 1,C < 0. 

V ' (1 - g4)a/4 ' ' - 



We show that y n = a n q n now satisfies 



q n - Q (-cy n y n+1 + q a )(-cy n yn-i + q Q ) = { [f_ y y r ^f_ cy 



(q a - y n )(q a - cy n )q a n even 

n odd. 



In Section 5 we relate these equations to a-q-Py. In Section 6 we show a stable method of 
computing the recurrence coefficients. 

2 Preliminaries 

We consider orthogonal polynomials on the exponential lattice 

L = {±q n \nEN}, 0<q<l. 

The orthonormality condition with respect to a weight wonL is, 

i 

p n (x)p m (x)w(x)d q x = S mn . (2.1) 

1 

The g-integral is defined as the sum 

/\ oo oo 

f(x)d qX = (i - q ) f(<i k M k + (i - 9) E f{-q k )q k - 

1 fc=0 k=0 

The g-difference operator D q will play an important role in our results: 

D q f{x)={ x(q-l) ltX ^ U 
(f'(0) ifx = 0. 

We will consider even weights w(x) = w(—x). This implies that the orthonormal polynomials 
p n associated with the weight w will satisfy the symmetry property 

p n (-x) = (-l) n p n (x). (2.2) 
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The orthonormal polynomials of even degree will be even, those of odd degree will be odd. The 
recurrence relation will then take the form 

xpn(x) = a n+1 p n+ i(x) + a„p„_i(x), (2.3) 

i.e., all coefficients b n are equal to zero. 

We will use the following technique (see, e.g., [TT] ) to identify g-discrete Painleve equations: 

• Given an even weight w, let p n denote the orthonormal polynomials with respect to w. 

• Find the Fourier expansion of the polynomial D q p n . 

• Compare the coefficients of this polynomial and its Fourier expansion. This will give rise 
to a set of equations. 

• After a change of variable, we find a g-discrete Painleve equation. We check the property 
of singularity confinement, and try to find out what happens when q — > 1. 

3 Modified g-Freud polynomials 

We consider the weight 

on the g-exponential lattice L, with a > — 1. The g-Pochhammer symbol is defined as 

n— 1 oo 

(«; q)n = JJ (i - aq 3 ) , (a; q)oo = JJ (1 - aq 3 ) . 

j=0 j=0 

Observe that, in terms of the g-exponential function E q (z) = (—z; g)oo, w e have 

w (t/l^x} = \x\ a E q , (-(1 - g 4 )gV) 



and hence w(y/l — q A x) — > \x\ a e~ x when g —> 1 (see, e.g., [6]). This limit relation is the only 
reason for the presence of the (constant) denominator in w. So this weight can be called a 
g-analogue of the modified Freud weig ht \x\ a e~ x . It is easy to check that this weight satisfies 
the Pearson equation 

f x\ (1 — x A ) , . .„ „. 

w \-j = { ~^w(x). (3.2) 

Let p n be the orthonormal polynomials associated to w as in (12. ip . Because of (12. 2p we can 
put p n (x) = 7 n x ?1 + S n x n ~ 2 + and we have the following 

Lemma 3.1. For orthogonal polynomials with an even weight one has 

n—X r 
Tn-1 , 2 d n 

a n = ana — > a,- = — . 

In ~[ In 

Proof. These two equations can be obtained by comparing, respectively, the x n+l and x n ~ 1 
terms in the recurrence relation (12. 3p . □ 
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Lemma 3.2. For even n, the polynomials p n with weight 113.1]) satisfy the following structure 
relation: 

B A 

DqPn(x) = —-?-p n _ 1 (x) + —^Pn-zi.x) (3.3) 

I — q 1 — q 

with 

A n = a " a "- ia r 2 (3.4) 

qa+n— 3 v ' 

(n+1 n-2 \ 

E-J-^E-J ■ < 3 - 5 > 
3=1 3=1 / 



with 



For odd n, the polynomials p n satisfy 

B A 

D q p n (x) = — — p n -i(x) H — pn-zix) + lower order terms (3.6) 

1 — q 1 — q 

^ = a n a ia 2 _ 

/n+l n-2 \ _ a 

«-^(EsWX>?) + ^. (3.8) 

Remark 3.3. Notice the difference between the closed expression of the structure relation for 
even n, and the presence of all even lower degree polynomials for odd n. Despite this difference, 
the resulting Painleve equations for even and odd n will have a very similar structure. 

Proof. If we expand D q p n into a Fourier series, we obtain 

n-l 

D q p n {x) = ^2a jin pj(x) 

3=0 

with 

»i 



a hn = J D q p n (x)pj(x)w(x)d q (x). 
The symmetry relation (12.21) shows that aj >n = if n — j is even. For n — j odd we get 

oo 

2(1 -q)Y l D q p n (q k )p j (q k ) W (q k )q k 

k=0 

oo 

-2 £ ( Pn {q k+l ) - p n {q k )) PjtfMq") 

oo oo 



k=0 k=0 



Since n — j is odd, (12. 2p implies that both sums are finite. Now we perform a shift in the 
summation index of the first sum, so that it contains the expression w(q k ~ 1 ), and we apply the 
Pearson equation (13.21) on this. We can recognize all sums as g-integrals, and we obtain 



q " f w(x) q ° /J 



a>j,n = ~y / Pn(x)p j (x/q) :: ^— L d q x + ^ / p n (x)pj(x/q)x 3 w(x)d q x 

i ~ 1 J-i x i ~ Q J-i 

1 f 1 w(x) 

+ ~ / Pn{x)p 3 \x) d q X. (3.9) 

1 — q J_ 1 x 



From now on we have to make a distinction based on the parity of n. 

Case 1: n is even. Since n — j is odd, we know from (I2.2p that Pj(x) and pj(x/q) are odd 
polynomials, hence and Pj are polynomials, and the orthogonality relations (12.11) imply 
that the first and the third integral in (13.91) vanish. To obtain the value of the second integral, 
we write pj(x/q)x 3 as a linear combination of the orthonormal polynomials Pk{x): an easy 
calculation shows that 

Pi(x/q)x 3 = — — — :p j+3 (x) + ( — — ^ J+3 — r | pj + i(x) + lower order terms 

Using Lemma [3. II and the ort honor mality relations (12. ip we obtain (I3.3l) -( l3~31) . 
Case 2: n is odd. The second integral in (13.91) can be computed in exactly the same way as in 
the previous case. However, this time the contributions of the first and the third integral do 
not vanish. For the third integral, we need to write p n (x)/x as a linear combination of the 
orthonormal polynomials Pk{x). This can be done because n is odd, and hence p n is an odd 
polynomial. This yields 

1 '' ' ' n -p n -i(x) + ( JrAi i | p ra _ 3 ( x ) _|_ l ower order terms 



X 7n-l \7n-3 7n-l7n-3, 

As for the first integral in fl 3 . 9 j) . we can use the orthogonality relations (12.11) to write 

1 t \ / i \ w ( x )j f 1 , \Pj( x /v) -pj( x ) +Pj( x ) , u 

Pn{x)Pj{x/q) d q x = / Pn{x)— w(x)d q x 

1 x J _i X 

1 ( \ I \ W ^ X ) A 

p n {x)p j {x) d q x, 

1 x 

which is the same as the third integral in (13.91) . Combining these results with Lemma 13.11 we 
obtain <^M>- §£E§. □ 

Now we can use these structure relations (I3.3p -( l3~8l) to obtain relations between the recurrence 
coefficients a n . Comparing coefficients of x n_1 and x n ~ 3 in these structure relations and using 
Lemma 13.11 gives 



a n a n-l a n-2 ~ 1 + 



n-2 

^(l-g^aj-fl-fX,] (3.10) 



and 



for even n, and 



a n a n+l 



n-2 

2 



+ al + aU + (1 - g 2 ) ^ a) ) = (1 - <f )<f (3.11) 



n-2 



-a-n+3 2 2 2 

u n u n-l a n-2 



q n -\l - q 2 )^] ~ (?*° " <f~ 2 ) <-i (3-12) 



and 

n-2 



< \al+ x + al + aU + (1 - q 2 ) g a 2 j = (q'* - q n ) q^ (3.13) 

for odd n. The aim is to obtain a Painleve-type recurrence relation between a n+ i, a n and 
a n _i. Hence to get rid of the a n _2 in (13.101) and (I3.12p . we replace n by n + 1. Keep in mind 



that this changes the parity of n. The sum X^=i a j that ar i ses m this way, should then be 
seen as Y^j=i a j + a n-i> "where we can use (13.131) and ( 13. lip respectively to write these sums 
as a function of a n+ i, a n and a n _i only. These manipulations yield the following recurrence 
relations: 

< K+i + q- a ~ n+1 a 2 n + + <T 2n - a+ X+i«X-i) = (1 ~ <? n ) <? Q+n " 1 , " even (314) 



a n ( a ra+l 



+ <T n+ X + + g~ 2n " a+3 an+ian«n-i) = (<T" ~ ? n ) q a+n ~\ n odd. 



This is a g-deformation of the discrete Painleve I equation d-Pi, which in its most general form 
is given by 

an + b + c(-l) n 

x n+ i + x n + x n -i = ha. (3.15) 



Indeed, if we take x n = o? n j y/l — q A and we let q tend to 1, we get 

2n + a -a(-l) n 

oX n 

Putting y n = a^q 1 '™, we can factorize the equations: we obtain 

a n-<* {v v 1 __|_ n ^ (v v .^.f?"- q~ a yl for even n , , 

q [VnVn+l + q ) {VnVn-l +9 ) - | \ - y 2 for odd 71. 

If we write u n = q~ a y2n and v n = —y2 n +i, substituting n by 2m resp. 2m + 1 in ( 13. 16ft . we get 

f q 2m (1 - u m v m ) (1 - u m v m -i) = 1 - u 2 m , s 

U 2m+1+a (l-^ m )(l- Mm+1 t; m ) = ^" A ' J 

This set of equations could therefore be called a set of asymmetric g-discrete Painleve I equations 
(a-q-Pi). However, in the next section we will obtain a more general form of it. Concerning 
the asymptotic behaviour of a n as n tends to infinity (with fixed q), it follows easily from ( 13. 14ft 
that a n tends to 0, and that 



a 2 a 2 
lim -r^-r = lim y 2n = q a and lim = lim T/ 2 n+i = 1, (3.18) 



or, equivalently, 



lim u n = 1 and lim v n = —1. 



The set of equations (13. 16ft has the singularity confinement property: if y n = 0, then a singu- 
larity occurs for y n+ i but this singularity does not propagate infinitely: e.g. for n even, putting 
y n = e we obtain 



Vn 


= e 




Vn+l 


= -q a \ 


;i - q~ n : 


Vn+2 


= q a - 1 


(l-q~ n 






- q~ a ~ 2 



1 

— y n -iq 
e 

I + y^zi 

e q 

2\ 



Vn+3 = * : \ e + 0(e') 
1 — q n 



1 — n n 

Vn + * = « a+ \_ q nl + s yn-l + 0(e) 
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so the singularity is confined to y n+ i, y n +2, 2/n+3j furthermore y n+4 depends on the value y n -i 
before the singularity. The same holds for odd n: 

Vn = £ 

Vn+1 = {q- n -q a )- £ -y n -iq- n - a + 0{e) 

y n+ 2 = ( q a - 1 - q - n - 1 ) 1 - + y -^ + o(s) 

e q 

Vn + s = q a f^e + 0(e-) 

There are two (one for each parity) more critical cases that might give rise to singularities, 
namely when y n y n -i + q a = 0. Running the singularity analysis gives, for even n, 

q a 

y n = h e 

Vn-l 

y n +i = 5 1 h °y £ ) 

Vn-l e Vn-l 



The coefficient ofl/e however, is 0(e) itself as y n -\ satisfies 

q n ~ a ~ l (y n y n -i + q a ) {yn-iVn-x + q a ) = i - yLi 

or, after substituing y n , 

1 - 2/n-i = syn-\q n ~ X ~ a {y n ~iyn-2 + q a )- 
So there is no singularity to confine. 

One could hope to use the recurrence relations (I3.16P to compute the recurrence coefficients. 
Given y = 0, it would be nice to show that there is a unique solution to ( 13.161) which is positive 
for all n > 0. The existence of such a solution is clear since y n = a^ l /q n ~ 1 satisfies the given 
recurrence relations and is obviously positive for n > 0. Its uniqueness however can only be 
expected based on numerical experiments, but we have no proof of it yet. The method used 
in Section [6] is probably the key to proving this uniqueness. Using (11. ip and the g-binomial 
theorem, it is clear that 

= a2= /V^gWUxIHa: = (g^ 1 ;^ 

Vl ai ~ /i (^V; g 4 )oo \x\°d q x ' (q a+3 ; q 4 )^ ' 

The knowledge of y and y\ and the recurrence relation (13.161) allow us to compute the y n 
recursively. However, this turns out to be an unstable method. We computed y n to a certain 
accuracy for some particular choices of q and a. In Figured] we plotted the values of log|y n |, 
calculated with an accuracy of 200 digits, with parameters q = 0.9, a = 5. For this choice of 
parameters, the computed values seem to satisfy the limit behaviour as stated in (13.181) up to 
n ~ 90, but for larger n the sensibility of the non-linear equations on the initial values destroys 
this behaviour. This is why Section will be devoted to a more stable way of computing the 
recurrence coefficients. 
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Figure 1: The result of computing log \y n \ from (I3.16P with a = 5, q = 0.9, with an accuracy 
of 200 digits. The odd-even dependence is clearly visible. 

4 Another modified discrete g-Freud case 

In this section we consider a generalization of the previous case: let the weight be given by 

, s = \A a {(fx 2 \q 2 ) 00 {cq 2 x 2 -q 2 ) 00 

y ] (1 - g 4 r /4 1 ' } 

with c < 0. This weight satisfies the Pearson equation 

(x\ (\-x 2 ){\-cx 2 ) , . 
W {q) = f W(X) - 

If 

c= -l + a^l-g 4 (4.2) 

then 

w (^Jl -g 4 x) -> \x\ a e- x4 - 2ax2 as q -> 1, 

so this weight is a g-analogue of the Freud weight \x\ a e~ x4 ~ 2ax2 . Obviously the choice c = — 1 
gives us the case considered in the previous section. The choice c = leaves us with a g-analogue 
of the weight |x| Q e _3; , which was already studied in [TTJ in the case a = 0. 



Lemma 4.1. For even n, the orhogonal polynomials p n for the weight w in (f^.ip satisfy the 
following structure relation: 

B. A 

D q p n (x) = — —Pn-^x) + — —p n ^ 3 (x) 

1 — q 1 — q 

with 

An ~ ~° q a+n-3 
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n+l n-2 

2 



' V j=i j=i 



For odd n, the polynomials p n satisfy 



with 



B A 
D q p n (x) = - — — p n -i(x) + - — — p„_ 3 (x) + lower order terms 



/i — a\ a n-l 




■^ + (c+l) ° 



qa+n—1 



Proof. The proof follows the same steps as the proof of Lemma 13.21 The only difference is 
that in the Pearson equation, the quartic polynomial 1 — x 4 has now got to be replaced by 
1 - (c + l)x 2 + cx 4 . □ 

Performing the same manipulations as in the previous section leads us to a recurrence relation 
satisfied by the a n : 



"1 



(c+l-c (a 2 n+1 + q- a - n+1 a 2 n + q 2 a 2 n _ x - cg-^+X+iaX-i)) = (1 ~ <f) q a+n ~\ n even 
a 2 n (c + 1 - c (a 2 n+1 + q- n+l a 2 n + q 2 a 2 n _ l - aT an - a+ X+i<M-i)) = {q' a - q n ) q a+n ~\ n odd. 

(4.3) 

Putting again y n = a^q 1 n we obtain a g-discrete Painleve I equation, of which (13 . 16[) is a 
special case: 

n-n i n\ i n \ \ (<7 Q — Vn) {q° — cy n ) q~ a , for even n ,. .v 

ft (-cy n y n+ i + <f) (-cy n y n ^ + g «) = ( [J _ c J"> for odd n _ (4.4) 

We can rewrite these equations into a more familiar form for asymmetric discrete Painleve 
equations: denoting u n = y2 n q~ a and v n = q/2n+i we get 

q 2n (1 - U n V n ) (1 - UnVn-x) = (1 - U n ) (1 - CU n ) , * 

q 2n+1+a (1 - u n v n ) (1 - u n+1 v n ) = (1 - v ft ) (1 - v n /c) • 

Besides q, this set of equations contains the parameters a and c, and it is the most general form 
of the asymmetric g-discrete Painleve I equations we obtain. Hence we call this set a-q-Pi. It 
is interesting to notice that these modifications w.r.t. the previous section 'survive' in the limit 

2 

case q — > 1: if we put x n = ^ — and we let q tend to 1, we obtain 



2n + a - a(-l) n 



j%r, 



which is still an instance of the discrete Painleve I equation (13.151) . Here a and c are related as 
in (14.21) . The introduction of the additional parameter c has no influence on the limit behaviour 
of the recurrence coefficients a n : ( 13.181) still holds. To show this we distinguish the following 
cases. We only mention the proof for even n, the odd n case being analogous. 

• c = — 1 is the case from the previous section. 
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> c > —1: for even n, (14. 3p gives 



4 < g 2a (l-g n ) 



g2n-2 _ c 



so a„ — > and ?/„ = is bounded. Denoting A = limsup?/™, (14.31) gives the quadratic 

equation A + 1 — = <7 Q . This equation has a negative solution A = q a /c which 

we can reject, and a positive solution A = q a . The same argument can be used on 
B = liminf y n , and one obtains (13.181) . 

c < —1, n even: from (14.31) we get 

^ 2 / a? 



. . c+l-c— ^- < l-q n < 1. 

qCt+n— 1 \ qa+n—1 I — ' 

a 2 

Looking at this as a quadratic inequality for y n = we obtain 



and hence y n G (1/c, 1). A similar argument with limsup and liminf yields (13.181) . 

• c = 0: then the result follows immediately from (14.31) . 

Once more, the singularity confinement property is fulfilled: taking y n = e for an even n, one 
finds that y n+ i and y n +2 ar e 0(l/e), y n+ 3 is 0(e) and 

= ^^^1^ [q/n-l(l " <f ) - (1 + C)(l - g^)] + 0(6). 



Now it is not obvious that the constant term is nonzero: if y n -i = ^cn-^) — ' then this is not 
yet a proof that the singularity confinement property holds. For odd n a similar argument 
yields 

Z/n + 4 = c(1 ^"J n+3) [cy n -i(l ~ q n+a ) - (1 + c)(l - <T V] + 0(e) 

and a problem could arise if y n -\ = 9 ^(i-g^+^) " • However (consider the case n is even, the odd 
n case being completely analogous) the zero at y n+4 gives rise to singularities in y n+ s, y n +6, Vn+7, 
but the singularity vanishes in y n+ s- Moreover, for the same reasons as above, this y n+8 can 
only be zero if y n+ 3 = fEi^R (which is nonzero since \q\ < 1 and c ^ -1, the case c = — 1 
being considered in the previous section), while the computation gave y n +z = 0(e). So this 
assures us that even in the worst case, the singularity is confined to y n +i, ■ ■ ■ ,Vn+7- Again, as 
in Sectional there are no singularities arising from —cy n Vn-i + q a = due to fine cancellations. 
Concerning the use of the recurrence relations to compute the recurrence coefficients a n of the 
orthogonal polynomials, we now start with y$ = and 

EOO qk(a+3) 
k ~ 



_ 2 _ ^fc=Q {q 2 ;q 2 )k{cq 2 ,q 2 )k 

yi — «i — ^ gMa+i) • 

Z^fc=0 {g 2 ]q 2 )k{cq 2 tq2)k 

The same remark as in the previous section holds: this method is unstable, as is shown in 
Figure El We refer to Section [6] for a stable computation method. 
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10 20 30 40 



Figure 2: The result of computing log \y n \ from 04. 4p with a = 2, q = 0.5, c = —1/3, with an 
accuracy of 200 digits. 

Remark 4.2. We restricted ourselves to the choice c < 0, with the case c = being a limiting 
case where the Painleve equations give explicit expressions for e.g. y n . However, numerical 
experiments strongly suggest that also for < c < 1 (in which case the weight w(x/q) is still 
positive on [—1,1]), the results on the limit behaviour A3.18\) still hold, and the singularity 
confinement property is still valid as well. Even for c > 1, where for small n we observe that y n 
can be negative, we still observe that y n is bounded and oscillating, however the even and odd 
subsequences do not tend to q a and 1 respectively anymore, but to q a /c and 1/c, respectively. 

5 Connection to a-q-Py 

It is worth noticing that the obtained recurrence relations can be seen as a limiting case of the 
set of asymmetric discrete Painleve equations a-q-Py- This set can be found as the discrete 
Painleve equation connected to the affine Weyl group E\ in Sakai's classification(see, e.g., jl]), 
or as a-q-Py in [llj. The equations are 



,.. W1 % (u n - l/p) (u n - 1/r) (u n - 1/s) (u n - l/t) 
1 - u n v n ) (1 - U n V n ^i) = 



(1 - u n v n ) (1 - U n+1 V n ) 



{u n - bp n ) (u n - p n /b) 

On - P) On - On ~ «) On ~ *) 



On - aw n ) (v n - w n /a) 
with prst = 1. Consider this system with the particular choice of parameters 

p = 1, r = c, s = k, t = — , b = ck, p n = q 2n , a = K,w n = q 2n+a+ 

en 

Then, letting k tend to 0, we obtain as a limit the set of equations ( 14.51) . 
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6 A stable method for computing the recurrence coeffi- 
cients 

A more stable way for computing the recurrence values is by writing (I4.4p as a system of 
quadratic equations in u n := y% n and v n := ?/2n+i, in the same way as we obtained (14.51) : 

u 2 n q- a (c 2 q 2n v n v n ^-c)+u n (-cq 2n (v n + 6 n _i) + c + 1) + q a (q 2n - 1) = 
v 2 (c 2 q 2n+1 - a u n u n+1 -c)+v n (-cq 2n+1 (u n + u n+1 ) + c+l) + (q 2M = 0. 

Computing the discriminants and opting for their positive roots (as we know, the u n and v n 
are related to the recurrence coefficients which are positive) in the expression of u n and v n , we 
find 

u n = fn(-cq 2n (v n + t> n _i) + c + 1, c 2 q 2n v n v n -i - c) 
v n = 9n(-cq 2n+1 (u n + Un+i) +c+l, c 2 q 2n+1 - a u n u n +i - c) 
with 

, s -x +v / x z+4(l- q ^)y 
Jn\^jy) — 2q~ a y 

, N -x+ x /x 2 +Ml-q 2n + 1 + a )y 

9n{x,y) = y Ty . 

We now define an operator T, acting on the space of double non-negative rows: 



T 



where 
and 



t = (£o,6,..o \ = re 

,2n/„ 1 „ \ 1 1 1 „2„2n 



£0 = °, i' n = f n (-cq n (r]n + Vn-i) + c+l,c q n rjnVn-i - c) for n > 



Vn = 9n(-Cq 2n+1 ^ n + ^n+l) +C+1, C^^Un+l ~ c) for U > 0. 

The solution of the equation which comes from our recurrence coefficients, denoted by (u,v), 
now coincides with a fixed point of this operator T. We know that (u, v) > (0, 0) (inequal- 
ities should be interpreted termwise: as holding between any two elements on corresponding 
positions). Unleashing the operator T (on any inequality we will use) reverses inequalities. 

Lemma 6.1. If < (a,b), thenT(£,rj) >T(a,b). 

Proof. It is not hard to check that all partial derivatives of /„ and g n are negative in the region 
where we need them. Which region this is, depends on the parameter c. □ 

We now construct the following sequence by repeatedly applying T to (0, 0). This sequence has 
increasing subsequence T 2fc (0, 0) and decreasing subsequence T 2fc+1 (0, 0). We denote their limits 
by (£,~,V~) an d (£ + , ? 7 + ), respectively. By continuity of T we get that T(£~,r]~) = (£ + ,?7 + ) 
and T(£ + ,77 + ) = (£ _ ,77~). Furthermore, any fixed point (£*,?7*) with positive components 
has to obey (£~,?7~) < — If one is a ble to show that (£ _ ,?7~) = ((, + ,t] + )> 

one proves there is only one fixed point, which necessarily coincides with the solution of the 
discrete Painleve equation coming from the recurrence coefficients. We did not find this proof 
(yet), but numerical experiments strongly suggest that indeed (£~,??~) = (£ + ,^ + )- This would 
also be enough to prove the uniqueness of positive solutions to (14.41) or (13. 16[) with y = 0, as 
mentioned in Section [31 See Figure [3] for the result using a particular choice of parameters, and 
Figure H] for a comparison between the two methods. 

Remark 6.2. If c = then &4-4\ ) immediately gives the exact recurrence coefficients: y n = 
q a — q n+2a for even n, and y n = 1 — q n+a for odd n, which gives 

g n+a-1 (l — g n+a ) for even n 



a . 



,2 

n ^1 „n — 1 1 



1 — g n+Q ) for odd n. 



13 



0,39- o 



0,38- 



0,37- 



0,36- 



0,35- 



o o o o o 



n — i — i — i — r 
2 



n — r - 
4 



10 



Figure 3: Approximation of the recurrence coefficient y\ for a = 2,q = 9/10, c = —1/2; on the 
horizontal axis is the number of times we applied the operator T to the starting value (0, 0). 



jp H SB SB E3 



a i 



i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i 
2 4 6 8 10 12 14 



Figure 4: log \y n \ with a = 2, q — 1/2, c = —5/2, computed using the recurrence relations f)4.4p 
with 20 digits accuracy (circles), and the approximation obtained after applying once (squares), 
or three times (crosses) the operator T of Section [H] to the starting value (0, 0). 



14 



7 Conclusions 



We used a generalized g-Freud weight in order to find a g-discrete Painleve equation. What 
we found was an asymmetric form of d-Pi which has (as we believe) never been described 
before. The asymmetric form emerges in a very natural way from the parity of the orthogonal 
polynomials associated to the weight. We showed its relation to a — qPy and gave a stable 
method to compute the recurrence coefficients for the orthonormal polynomials associated to 
this weight. 

The technique used in this paper is not sufficient to handle weights that are not even. In these 
cases, one has a non-zero b n entering the system, which could lead to new results. This is a 
possible interest for future research. 
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